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Abstract 

Photometric observations of the eccentric eclipsing binary VI 143 Cygni were performed 
during the August-September 2000 and July 2002, in B and V bands of the Johnson system. 
The analysis on both light curves was done separately using the 1998 version of Wilson's 
LC code. In order to find a new observed rate of apsidal motion, we followed the procedure 
described by Guinan and Maloney (1985). A new observed rate of apsidal motion of 3.72 
degrees/lOOyr was computed which is close to the one reported earlier by Khaliullin (1983), 
Gimenez and Margrave (1985), and Burns et al.(1996). 

1 Introduction 

VI 143 Cygni (HD 185912; HR 7484; BD+54°2193; V max = +5.86; B-V=+0.46; a = 
19 38 m 41 s .18; 5 = +54°58'25".7) is a double-lined eclipsing binary, consists of a pair 
of F5V stars with high orbital eccentricity (e=0.540) and a relatively long period of 
7.640 days (see Fig.l). According to Andersen et al. (1987) who have determined the 
stellar and orbital properties of this system accurately, the corresponding radii and 
masses are n = 1.346 ± O.O23i? , r 2 = 1.323 ± 0.023^ and m x = 1.391 ± O.O16M , 
m 2 = 1.347 ± 0.013M©, respectively. 

It is known that VI 143 Cyg is one of the best examples of eclipsing binaries with 
apsidal motion, in which the observed rate of apsidal motion is grater than the value 
predicted by general relativity and stellar evolutionary models. 



1 



Mini 




To Earth 



Figure 1: Relative orbit of VI 143 Cygni drawn in scale. 

The observed rate of apsidal motion is due to the contribution of two terms; a 
classical term as well as the general-relativistic term. In this sence, the observational 
apsidal motion rate is 

^obs = &cl + U GR , (1) 

where Co c i denotes the classical or Newtonian term and ujgr is the relativistic contribu- 
tion which can be determined using the formulas given by Gimenez (1985). In the case 
of V1143 Cyg, the observed rate of apsidal motion is u obs = 3°.52/10(F r ± 0°. 72/10(F r 
(Burns at al. 1996) while Andersen at al. (1987) calculated a faster theoretical ap- 
sidal motion of iotheo — 4°. 25/100^ ± 0°.72/100 ?/r in which the expected relativistic 
and classical (Newtonian) contributions to apsidal motion are uj c \ = 2°.39/100 yr and 
^gr = l°.86/100 yr , respectively. It is seen that the classical and relativistic contribu- 
tions are of the same order in this system. 

2 Observations 

VI 143 Cygni was observed during 24 nights from July to September 2000 at Biruni 
Observatory of Shiraz University (Longitude: 52°31' E, Latitude: 29°36' N). Obser- 
vations were made with a 51cm cassegrainian telescope equipped with an uncooled 
RCA4509 multiplier phototube. Two stars HD 184240 (V max = +6.30) and HD 186239 
(Vmax = +7.10) were selected as comparison and check stars respectively. The integra- 
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tion time for all of the observations were fixed to 10 seconds. The output signals of 
the photomultiplier were fed to a computer after amplification, using an A/D conver- 
tor. The measurements were made using B and V filters of intermediate-bandpass blue 
^■max = 4400A and yellow \ max = 5530A which are matched closely to the Johnson's 
UBV system. Times were converted to Heliocentric Julian Day Number (HJD). Data 
reduction and atmospheric corrections were done to obtain the complete light-curves 
in two filters, using a computer code developed by G. P. McCook. Figures 2 and 3, 
represent the observed light curves in B and V filters respectively. In order to enhance 
the accuracy of the present work , two other minima were observed on July 16 and 18, 
2002. This time, the star HD 185978 (F8; m v = +7 m 8) was selected as the comparison 
star. 
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Figure 2: Observed light curve of V1143 Cygni in the B filter. 



3 Times of minima and light curve analysis 

From the observed light curves, heliocentric times of minima (one primary and one 
secondary) were computed by fitting a Lorentzian function to the observed minima 
data points (see Figures 4 and 5). This function can be expressed as 

2A w 
n 4(x — x c ) 2 + w 2 



y = y» + — iuz - v i -.2 > ( 2 ) 



where y Q is the baseline offset, A is total area under the curve from baseline, x c is the 
center of the minimum and w is full width of the minimum at half height. The minima 
were calculated according to the ephemeris given by Andersen et al. (1987) 

Min.I = HJD2U9234MU + 7^.64075217 x E, (3) 
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Figure 3: Observed light curve of VI 143 Cygni in the V filter. 

and are given in Tables 2 and 3. The corresponding errors in primary and secondary 
minima are: ±0.00066 and ±0.00220, respectively. Meanwhile, in each filter, the depths 
of minima are as follows: 

Filter B: Min.I: m 53 ± 0.02, Min.II: m 25 ± 0.02 
Filter V: Min.I: m 48 ± 0.02, Min.II: m 23 ± 0.02 




Figure 4: A sample Lorentzian fit to the primary minimum. 

The B and V light curves of VI 143 Cyg have been analyzed separately by using 
the Wilson code in order to derive photometric elements of this system. The program 



4 



n 1 1 1 1 1 1 1 1 1 ■ 1 1 r 




PHASE 

Figure 5: A sample Lorentzian fit to the secondary minimum. 



consists of two main FORTRAN programs LC (for generating light and radial velocity 
curves) and DC (to perform differential corrections and parameter adjustment of the 
LC output). The model which upon the program is based on, has been described and 
quantified in papers by Wilson (1979, 1990, 1993). Since the system VI 143 Cygni is 
detached with both components residing well inside their respective Roche lobes (Burns 
et al. 1996), the solution was performed in mode 2. In our analysis, we assumed a 
value of zero for the third light (Z3 = 0). Also we fixed the ratio of the axial rotation 
rate to the mean orbital rate for stars 1 and 2 (Fi,F 2 ). Before running the LC code, 
the position of the preastron was estimated from the apsidal motion study as discussed 
in section 4. In order to optimize the observed parameters given in Table 1, we wrote 
an auxiliary computer programm to compute the sum of squares of the residuals in 
both filters by using theLC output. Figures 6 and 7 show the sum of the squared 
residuals versus selected parameters for the two filters. For example, Figure 6(a) rep- 
resents the weighted sum of the squared residuals J2 ujr ' 2 (SSR) versus eccentricity (e) 
and inclination angle (i) for filter B. 

Figures 8 and 9, show the optimized theoretical light curves together with the 
observed light curves in B and V filters. The theoretical curves correspond to the 
optimized parameters given in Table 1. 
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Figure 6: The weighted sum of squared residuals ujr ' 2 (SSR) versus parameters using 
Wilson's LC code in B filter. 



Table 1. 

Optimized parameters of VI 143 Cygni. 



Parameter 


Filter B 


Filter V 


i 


87.3 ±0.1 


87.1 ±0.1 


e 


0.536 ±0.005 


0.539 ±0.005 




0°.855± 0.004 


0°.855 ± 0.004 




0.99 ±0.01 


0.98 ±0.01 




17.4 ±0.5 


18.0 ±0.5 


n 2 


19.6 ±0.5 


20.0 ±0.5 


U 


0.89 ±0.05 


0.82 ±0.05 



Finally, the observed times of minima on July 16 and 18, 2002 are calculated ac- 
cording to the ephemeris 
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Figure 7: The weighted sum of squared residuals ^f 2 (SSR) versus parameters using 
Wilson's LC code in V filter. 



Min.I = //JD2447087.5669 + 7 d .64075095 x E, (4) 
given by Burns et al. (1996) and are tabulated in Tables 2 and 3: 

4 Apsidal motion 

Due to the deep, narrow eclipses of VI 143 Cyg and its high eccentricity, the rate of 
apsidal motion can be determined exactly by analysis of primary and secondary eclipse 
timings. In their paper, Guinan and Maloney (1985) described the procedure that must 
be followed to determine the apsidal motion rate from the change in the displacement 
of the secondary minimum from the half point (0.5 phase) according to 

D= [(i 2 -fi) -0.5 x Period]. (5) 
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Figure 8: 
data. 



Theoretical light curve (LC output) in filter B. Points are the observational 



D in turn, is related to uj the longitude of preastron by the formula given by Sterne 
(1939): 



^ P r 1 / e cos(u; 
D = — tan" 1 — 

vr L Vl-e 



2^/2 



+ 



ecos w 



1 



e 2 sin 2 fu;) 



e 2 ) 1/2 ], 



(6) 



where P is the period and e is eccentricity. The observed photoelectric times of primary 
and secondary minima from 1969 to 2002 are collected in Tables 2 and 3. As you 
can see from Table 4, the slow decrease in D is due to the advance of the line of 
apsides of the orbit. By computing the slope of a line fitted to all of the secondary 
minima given in Table 3, we determined an observed rate of apsidal motion of u b s = 
3°. 72/100^ ±0.37/100f r (see Figure 10). 
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Figure 9: Theoretical light curve (LC output) in filter V. Points are the observational 
data. 



Table 2. 

The photoelectric times of primary minima for V1143 Cyg. The O-Cs are computed 



according to equation 7. 



H.JD. 


O-C(day) 


Epoch 


Reference 


(2400000. +) 








39339.616 


-0.22773 


-376 


Snowden and Koch (1969) 1 


39385.6881 


-0.00011 


-370 


Snowden and Koch (1969) 


40837.4313 


0.00017 


-180 


Battistini et al. (1973) 


41135.4208 


0.00033 


-141 


Battistini et al. (1973) 


42212.7651 


-0.00142 





Koch (1977) 


42617.727 


0.00061 


53 


Koch (1977) 


43305.3943 


0.00021 


143 


Guinan et al. (1987) 


45253.7858 


-0.00008 


398 


Gimenez and Margrave (1985) 


47087.5669 


0.00049 


638 


Burns at al. (1996) 


48019.73800 


-0.00016 


760 


Caton and Burns (1993) 


49234.6144 


-0.00336 


919 


Lacy and Fox (1994) 


51771.34410 


-0.00338 


1251 


Dariush et al. (2001) 


52474.29443 


-0.00225 


1343 


Dariush et al. (2003) 



1. This minimum is observed spectroscopically and because of its 
unusually high 0-C, it is not shown in the 0-C diagram. 

2. This minimum was published as 2439385.6831 by Hamme and Wilson 
(1984) which seems to be misprint. 
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Table 3. 



The photoelectric times of secondary minima for VI 143 Cyg, together with 
the computed values of D and u, using equation 6. The O-Cs are computed 

according to equation 8. 



H.JD. 


O-C(day) 


Epoch 


D 


UJ 


Reference 


(2400000. +) 












38932.932 


-0.00431 


-430 


1.8685 


47.97 


Snowden and Koch (1969) 


38978.7807 


0.00000 


-424 


1.8727 


47.83 


Snowden and Koch (1969) 


42615.75 


-0.01896 


52 


1.8439 


48.82 


Koch (1977) 1 


43066.575 


0.00286 


111 


1.8646 


48.11 


Koch (1977) 


44487.7482 


-0.00002 


297 


1.8579 


48.34 


Gimenez and Margrave (1985) 


47085.5910 


-0.00598 


637 


1.8449 


48.78 


Burns at al. (1996) 


51792.28645 


-0.00122 


1253 


1.8370 


49.05 


Dariush et al. (2001) 


52472.30445 


-0.00834 


1342 


1.8281 


49.35 


Dariush et al. (2003) 



l.Due to its unusually high 0-C, it is not included in determination of 
the observed rate of apsidal motion. 
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Figure 10: The displacement of the secondary minima given in Table 3, from the 0.5 
phase versus epoch. The observed rate of apsidal motion (u b s ) can be calculated for 
the slope of this curve. 

The residuals given in Tables 2 and 3 are calculated for each of the minima according 
to the ephemeris given by Gimenez and Margrave (1985) 

Min.I = HJD2U22 12. 76652 + 7 d .64075217 x E, (7) 
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and 



Min.II = #J J D2442218.45092 + r.64073165 x E. 



(8) 



The computed (O-C)s versus Julian day is plotted in Figure 11 for both primary and 
secondary minima. This diagram shows no changes in period within the 0.01 days. 
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Figure 11: The O-C diagram in days for all the primary and the secondary minima 
tabulated in Tables 2 and 3. 



5 Results and discussion 

Table 4 contains independent determinations of the observed apsidal motion of VI 143 
Cyg together with the corresponding period of apsidal revolution U. The relation 
between the apsidal motion period U and the observed rate of apsidal motion {to b s ) 
has a simple form 

U = °™ (9) 

^obs 

where u b s is expressed in degrees per cycle and P is the anomalistic period expressed 
in days. To determine a more accurate value for u obs we need more accurate timings of 
secondary minima. From Table 4, it seems that expanding our observational baseline, 
may decrease the discrepancy between u; f> s and Utheo- Figure 12 shows D as a function 
of uj using equation 6 for VI 143 Cygni and DI Herculis. In the case of DI Herculis 
(P=10.550 days, e=0.49) the discrepancy is much larger than that of V1143 Cygni. For 
this system the observed rate of apsidal motion is about one-third of the theoretical one 
which is due to the classical and relativistic contribution (Claret, 1998; Dariush and 
Riazi, 2003). The dots in Figure 12, represent our observational period for the VI 143 
Cyg and DI Herculis apsidal revolution. It is clear that up to now, our observational 
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baseline covers only a very small fraction of the period of apsidal revolution U. In 
the case of DI Herculis, the presence of a third body is a plausible explanation for 
the discrepancy between u obs and totheo, but until no further observational evidence has 
been reported this possibility. Our results for VI 143 Cyg supporting 3°.72±0.37/10Ch/r 
is slightly closer to the theoretical apsidal motion rate (4°. 25 ± 0.72/10Ch/r), computed 
with previous studies of this system 
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Figure 12: D as a function of u using equation 6 for VI 143 Cygni and DI Herculis. 

Table 4. 



Determined rate of apsidal motion for VI 143 Cygni 



c^/lOO"- 


±Error 




Source 


3.49 


±0.38 


10320 


Khaliullin (1983) 


3.36 


±0.19 


10710 


Gimenez and Margrave (1985) 


3.52 


±0.72 


10230 


Burns et al. (1996) 


3.72 


±0.37 


9680 


Present study 1 



l.This value is computed from all the secondary minima presented 
in Table 4 except 2442615.75 of Koch (1977). Including this 
minimum, it is changed to 3.45 degree/lOOyr +/- 0.65 degree/lOOyr . 

Acknowledgments 

We would like to thanks Mr. Mehdi Nazem for his help during the observations. 




12 



References: 



Andersen, J., Garcia, J. M., Gimenez, A., and Nordstrom, B., 1987, Astron.Astrophys., 
174, 107. 

Battistini, P., Bonifazi, A., and Guarnieri, A., 1973, IBVS 817. 

Burns, J. F., Guinan, E. F., and Marshall, J. J., 1996, IBVS 4363. 

Caton, D. B., and Burns, W. C, 1993, IBVS 3900. 

Claret, A. 1998, Astron.Astrophys., 330, 533. 

Claret, A., and Willems, B., 2002, Astron.Astrophys., 388, 518. 

Dariush, A., Afroozeh, A., Riazi, N., 2001, IBVS 5136. 

Dariush, A., and Riazi, N., 2003, Astrophys. Space. Sci., 283, 253. 

Dariush, A., Zabihinpoor, S. M., Bagheri, M. R., Jafarzadeh, Sh., Mosleh, M., and 

Riazi, N., 2003, IBVS 5456. 

Gimenez, A., 1985, Astrophys. J., 297, 405. 

Gimenez, A., and Margrave, T. E., 1985, Astron.J., 90(2), 358. 

Guinan, E.F. and Maloney, F.P., 1985, Astron.J., 90, 1519. 

Guinan, E. F., Najafi, I., Zamani-Noor, F., and Boyd, P. T., 1987, IBVS 3070. 

Hamme, W. V., and Wilson, R. E., 1984, Astron.Astrophys., 141, 1. 

Khaliullin, Kh. F., 1983, Astron. Tnsk, No 1262. 

Koch, R. H., 1977, Astron. J., 82, 653. 

Lacy, C.H. and Fox, G.W., 1994, IBVS 4009. 

Semeniuk, I., 1968, Acta. Astron., 18, 1. 

Snowden, M. S., and Koch, R. H., 1969, Astrophys. J., 156. 

Sterne, T.E., 1939(a), Mon. Not. R. Astron. Soc, 99, 451. 

Sterne, T.E., 1939(b), Mon. Not. R. Astron. Soc, 99, 662. 

Wilson, R.E., 1979, Astrophys. J., 234, 1054. 

Wilson, R.E., 1990, Astrophys. J., 356, 613. 

Wilson, R.E., 1993, in New Frontiers in Binary Star Research, ed. K. C. Leung and 
Nha, I.S., A.S.P. Conf. Ser., 38, 91. 



13 



